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1. Introduction 



During tlie last two decades, diffusion-limited chemical reactions have attracted 
considerable interest. In particular, the one- and two-species annihilation reactions 
A + A^C and A + B ^ C are known to exhibit anomalous kinetics in lower dimension. 
For the one-species case, the upper critical dimension — the dimension above which the 
rate equation ria(t) = —kn'^{t) (where na{t) is the concentration of A particles) gives 
qualitatively correct results — is 2 [||, ^, ||, ||, |^. For the two-species case, the situation 
is more complex, because of the presence of a conserved quantity, namely the difference 
concentration of A and B particles. Depending on the initial conditions, the upper- 
critical dimension is 4 (for homogeneous conditions with the same initial concentration 
of both particle types) or 2 (for segregated initial conditions or non-equal initial A and 
B concentration, etc.) [jl], H, 0, ||, ||. Among other methods, those results have been 
derived using renormalization group techniques, drawing a rather complete picture of 
the different universality classes involved in these reactions ^, ||, 

In this paper, the question we want to address is the following: what happens if the 
backward decombination reaction is allowed with a given probability? In most physically 
interesting systems it is unlikely that this possibility is totally forbidden. Although for 
extremely small backward probability, one expects the effects to be very small (and 
even unnoticeable for not too large observation time), a fundamental change occurs in 
the system. Instead of decreasing toward a non-equilibrium steady-state, the system 
should eventually reach an equilibrium state. Moreover a new conserved quantity can 
be constructed reflecting the conservation of mass or energy: na{t) + 2nc{t) for the one- 
species reversible reaction A + A ^ C and Uait) + nh{t) + 2nc(t) for the two-species 
reaction A + B ^ C (where nb{t) and nc{t) are respectively the B and C particles 
concentration). 

At the mean-field level, the rate equations give an exponential approach toward the 



equilibrium state. However it was soon recognized ||T0|, [TT], [1^, |T3l that the conserved 
quantity obeys a diffusion equation, and thus its initial fluctuations should decay with 
a power law: t"*^/^ (d is the dimension of the system). In the long time limit, this 
power will always overcome the exponentially fast decay of the rate equations, for 
any dimension d > 0. The upper critical dimension is thus infinity. At first sight 
surprising, this result has been confirmed in various ways. For the single-species reaction, 
Zeldovich and Ovchinnikov |11| obtained the approach to equilibrium in the low density 
limit of a field theory in three dimensions. This result was extended later |jl3[ for 



arbitrary dimensions, but in the framework of an uncontrolled approximation. In the 
A + B ^ C case, the power law decay was first obtained in three dimensions [TDI from 
heuristic considerations. Later Kang and Redner |]12|, using an argument based on the 
fluctuations, and Burlastky et al [O, assuming a closure of the hierarchy, extended 
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this result to arbitrary dimensions, but with different amphtudes. Our purpose in this 
paper is to derive, using field theory techniques and Langevin equations, the asymptotic 
approach to equilibrium in a controlled and exact way. Whereas we are able to show the 
universality of the power law exponent, the amplitude proves to be model-dependent. 

The paper is organized as follows. In section ^ we consider the single-species 
A + A C reaction. From the master equation describing our model, we map the 
problem to a set of two Langevin equations for the random variables a and c, with 
complex noise. The concentrations na{t) and nc{t) are then given by the average of the 
random variables a and c over the noise: Ua = (a) and ric = (c). Exploiting the fact 
that na{t) + 2nc{t) is conserved by the dynamics, i.e. that a + 2c obeys a noisy diffusion 
equation, we can write down the long time behavior of its two-point correlation: 

{{a + 2c)2) - (a + 2c)2 ~ {c^ - c^){MDt)-'"^ (1) 

(coo is the steady-state density of C particles, Cq the initial one, and D the diffusion 
constant of both A and C particles). We can then set up a controlled approximation 
scheme to obtain the approach to the equilibrium, based on the fact that the previous 
correlation goes to zero as t grows to infinity. We find that 

(c) - Coo ~ .3 (00 - c^WT^Dt)-"' (2) 

(A and are respectively the rate of the forward and backward reactions and a^o is the 
equilibrium density of A particles). The exponent of the power law is universal, as it 
comes from the conservation law, but the amplitude is model-dependent. It is however 
interesting to note that Burlatsky et al |jl3| found the same result (up to a factor of 
2) for a slightly different model, their results relying, nevertheless, on an uncontrolled 
approximation (the closure of the hierarchy). Note that the final approach to equilibrium 
is governed by the sign of cq — Coo- In fact, the density is expected to reach its equilibrium 
value monotonically. 

In section ^, we consider the reaction A + B ^ C . Using the same method (in 
this case there is a second quantity for which the Langevin equation can also be solved 
exactly), we find 

(c) - Coo ~ ^[4? + /i' - A2(ao - 60)'] (co - Coo)(87rDt)-'^/2 

where ctab = A(ao + &o + 2co — 2coo) + /i, and Qq and bo are the initial A and B particle 
densities. An equivalent expression were obtained by Burlatsky et al ||T^ when = 6o- 



Similar conclusions concerning universality and monotonicity can be drawn. 

In section ^, various generalizations are considered. First we investigate the case of 
pair correlations between the reactants (i.e. A- A pairs for A + A ^ C and A-B pairs for 
A + B ^ C). We show that, depending on the initial fraction of correlated particles, the 
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approach to equilibrium can become non-monotonic, or even, in some special cases, the 
power law may change to a faster decay. The second part of this section is devoted to 
the interesting problem of unequal diffusion constants for the different species. Even in 
extreme cases (such as immobile C particles), the power law is unchanged and only the 
amplitude is modified. The last generalization we consider is the case of segregated initial 
conditions. When initially the A and B particles are spatially separated, a reaction front 
will develop. One natural question to ask is to which rapidity the width of this front will 
grow. We are able to confirm the scaling results obtained by Chopard et al [|T^] showing 
that this increase is governed by the mean-field behavior w{t) ~ t^^^. This result is also 
valid for immobile C particles. Final remarks are made in section ^. 



2. The A + A^C reaction 

2.1. The model and the formalism used 

Our starting point is the following continuous time master equation 
|p({m},M;t) 

= im[("^e, + l)P{- . . , mi - 1, + 1, . . . , {n}- 1) 

- miP({m}, {n};t)] 

+ ^ E EtK. + l)P({m}, . . . , - 1, ne, + 1, . . . ; t) 

-n,P({m},M;t)] 

+ Ao EtK + 2)(mi + . . , mi + 2, . . . , {n}; t) 

i 

- mi{mi - l)P({m}, {n}] t)] 

+ /i ^[(n, + l)P({m}, . . . , + 2, . . . ; t) - n,P({m}, {n}; t)]. (4) 

i 

The set {m} {{n}) denotes the occupation numbers of A {C) particles in each lattice site 
and P({m}, {n}; t) is the probability to find the configuration {m}, {n} at time t. This 
equation describes the evolution of the probability P in time. A given configuration can 
change due to four processes: by diffusion of A particles (second and third lines of (4), 
where Da is the diffusion constant of the A particles and £ is the lattice constant); by 
diffusion of C particles (with diffusion constant Dc, fourth and fifth lines). It will also 
change when two A particles merge into one C, with a microscopic reaction rate Aq (sixth 
and seventh lines) or, as shown in the last line of (4), when a C particle reacts producing 
two A (the corresponding rate is denoted by /i). In the diffusion terms, the sum over Cj 
stands for a summation over all the nearest neighbors of site i. In this respect, equation 
(4) models then the time continuous evolution of the reversible reaction A + A^ C on 
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a (i- dimensional hypercubic lattice, allowing for multiple occupancy on each site. 

For the time being, we choose the initial conditions to be given by an uncorrelated 
Poissonian distribution on each site and for each species: 



-rrii ~n. 
-ao-co TT Co 



P({m},M;0) = e----n^^' (5) 



where oq (cq) is the average occupation number per lattice site for the A (C) particles. 
A correlated initial condition will be considered later, in section ^. 

More than twenty years ago, Doi [|15| (see also [0) developed a procedure mapping 



the master equation to a second quantized representation by introducing sets of creation 
and annihilation operators. In turn this second quantized form can be mapped to a field 
theory (see [0). By now these various steps are well known and we shall only quote 



the results we need. 

Let us introduce the (complex) fields a, a, c and c. In terms of these fields, the A 
particle density is given by na(x, t) = ((a(x, t))) where ((■)) denotes an average over e~'^, 
where S is an action. In general, 

{{A[a,c])) = jvaVaVcVcA[a,c]e'^ I jvaVaVcVce~^ . (6) 

The script V denotes functional integration and S is the action corresponding to our 
reaction, obtained by the mapping of the master equation (4) 

S = j d'^xj^' dt [a{dt - DaV^)a + c{dt - DcV^)c + (a^ + 2a - c) (Aa^ - fic) 

- 6{t){aoa + cqc)], (7) 

where A = Aq^"^, cto = ^o/i'^ and cq = co/i'^. For convenience, the continuous space limit 
{i 0) has been taken. However an equivalent result can be written keeping the lattice 
structure. Note that any observable A can always be written as an expression which 
solely depends on the fields a and c and not on the response fields d and c (this comes 
from probability conservation). The double bracket notation in (^) stresses the fact 
that average is taken both over the dynamics (first line of (7)) and the initial conditions 
(second line). 

The density correlation function for the A particles is given by 

a(x,t) = (([a(x,t)+5'^(x)]a(0,t))). (8) 

Similar relations hold for the C particles. 

The analysis of the action using the renormalization group formalism has prowed 
to be extremely powerful for various types of reaction (see 0, ^ ||, Q for some recent 
examples). In particular it is especially well suited to distinguish universality classes. 
However, in our case, this analysis is not well adapted because the upper critical 
dimension is infinite and an expansion around it will clearly fail. In fact, to treat this 
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problem, we shall prefer to use the formalism of Langevin equations. These equations 
can be very easily obtained by replacing the quartic piece of the action with an integral 
over a noise variable: 

/+00 
dC cxp[aC - \(^/{^lc - \a^)] (9) 
-oo 

and integrating out the response fields a and c. One then obtains the two following 
equations 

{dt - DaV'^)a{x, t) = -2Aa(x, tf + 2/xc(x, t) + C(x, t) (10) 
(dt - L>cV')c(x, t) = Aa(x, t)' - //c(x, t), (11) 

where C is a complex Gaussian noise with zero mean value, whose correlation is given 

by 

(C(x, i)C(x', t')) = 2(/.c(x, t) - Aa(x, ty)S^^\^ - ^')S{t - t'). (12) 

Here, the single bracket notation stands for the average over the noise. There is no 
more need to average over the initial conditions as it has been explicitly done when 
integrating over the response field. 

Note that by using (11), one has, for homogeneous initial conditions, 

(C(x,t)C(x',0) = -2dMt))S^'\^-^')S{t-t'). (13) 

As the density is expected to reach a (reversible) stationary state, the noise correlation 
should vanish in the long time limit (lim^^oo dt{c{t)) = 0). This, together with the fact 
that the equation for the variable c(x, t) comes without explicit noise (except through 
the a dependence) are the central points of our analysis. 

These Langevin equations look very similar to the rate equations, with the addition 
of the noise. One then might question the necessity of deriving them in a such 
complicated way, as it would have been easier just to add noise to the rate equation. 
In fact this method allows us to derive exactly the noise and its noise-noise correlation 
function. In particular one sees that the noise is complex, a result which is at odds 
with the usual guess made when writing heuristic Langevin equations. It should also be 
emphasized that the variables a(x, t) and c(x, t) do not represent the density, because 
they are complex. The mean density na{t) {nc{t)) is given by the average of a(x, t) 
(c(x, i)) over the noise. One easily convinces oneself that this average give a real value 
for the density. 

An immediate consequence of the vanishing of the fluctuations at equilibrium is the 
possibility of computing the actual values of the equilibrium densities, which are simply 
given by their mean field values. Denoting by Uoo and Cqo such densities, one has 

Aa^ = I^Coo, (14) 
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which together with the conservation law 

ttoo + 2coo = ao + 2co, (15) 

give us 

aoo = + 8A(ao + 2co)//i - l). (16) 

It is easy to check that this resuh is a solution of the detailed balance condition of 
the master equation (4), which should hold as the stationary state is an equilibrium 
state. Note that we do not expect this result to be universal, i.e. apply for all models 
describing a reversible A + A ^ C reaction. In fact this result strongly relies on the 
multiple occupancy property and single site reactions of our model. It can be indeed 
shown that, in one dimension, a spin chain model of this reaction with exclusion 
process leads to a different steady-state density (equation (0) is replaced by another 
condition). However in the small density limit (dilute gas) this latter result converges 
toward the mean field value which is expected to be universal in this regime. 



2.2. Conservation law 

The next step in our analysis of the model is, of course, to obtain the approach toward 
equilibrium. It is easily seen that the rate equations give an exponential decay. The 
fluctuations are expected, however, to change this law. Our starting point to analyze this 
problem will be the two equations (10) and (11) which, with the noise-noise correlation, 
describe completely our model. In order to simplify, we shall consider now the case of 
equal diffusion constant Da = Dc = D. The case for which Da 7^ Dc will be considered 
in section ^ 

First one remarks that the quantity x = o + 2c obeys a noisy diffusion equation: 

(at-D,V2)x(x,t) = C(x,t) (17) 

This reflects the fact that the quantity na(x, t) + 2nc(x, t) = (x(x, t)) is conserved 
by the dynamics, which in turn is a statement about mass (or energy) conservation. 
Equation (|l^ is easily solved, and one finds 

X(x, t) = dt'J d^x' Go(x - x', t - t')C(x', t') + J dV G'o(x - x', t)x(x', 0) (18) 

where x(x, 0) is the initial condition and Go(x, t) is the free propagator: 

Go(x,t) = 9{t) (47rDt)-'^/2exp(-^) (19) 

{9{t) is the usual Heaviside step function). Using {Q = 0, 

(X(x,t)2) - (x(x,t))2 = -2 fdt, [8nD{t-t^)]-'^/'dMti)), (20) 

Jo 
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with (x(x, t)) = tto + 2cq. Although the exact structure of (c(t)) is not known (as this 
is precisely the quantity we want to compute), we do not need it to obtain (x(x, t)^) for 
long time, as, when t oo 

(X(x,t)2) - (x(x,t))2 = -2[8nDt]-''/' / dt, 5,(c(ti)) = -2{c^ - co)[87rDt]-''^\ (21) 

Jo 

At this stage several remarks have to be made. First, the integrand of expression 
(20) diverges when ti — > t. For d > 2 the integral is thus singular. This divergence 
is however artificial as it comes from the continuous space limit we took when writing 
the action, and it can be avoided by putting a short distance cut-off of the order of 
i in the space integration. In turn, once the integration over space is performed, this 
small distance cut-off will produce a cut-off function Ccfii"^ / D{t — ti)) which multiply the 
integrand of equation (|20|). The exact form of this cut-off function is unimportant. It 
should be a rapidly decreasing function for large x and it should go to 1 when x goes to 
(for example a possible candidate could be Ccf(x) ~ exp(— x^)). In the following, such 
regularization will always been understood when facing ultra-violet divergent integrals. 
A second remark concerns the sign of the variance of x which can be positive or negative, 
depending on the initial densities oq and cq, i.e. on the sign of cq — Coo- (Note that Coo 
depends both on qq and Cq and it is always possible to adjust Oq such that Cqo < Cq, or 
Coo > Cq.) The possibility of having negative variance plays a central role. In the next 
subsection, we shall show that the approach toward equilibrium is monotonic, at least 
for times at which our analysis apply. Whereas this result is not really surprising as it 
was already obtained by the rate equation, one should note that it refutes the results 
of a common method of solving such problems. From the two Langevin equations (10) 
and (11), a natural approximation for the density would be to consider the standard 
rate equations 

{dt - DV^)d{x, t) = -2Aa(x, tf + 2/ic(x, t) (22) 
{dt - L'V^)c(x, t) = Aa(x, tf - /ic(x, t) (23) 

with Poissonian random initial conditions. Whereas this approximation gives perfectly 
good results, both for short and long times (but not for intermediate times), for the 
two species annihilation reaction A + B ^ (/} (see for example |1^), in our case it 
fails to predict the anti-correlation of the conserved field. Indeed, it is easy to solve the 
equation for x = a + 2c. One has 

X(x, t) = I dV G'o(x - x', t)x(x', 0). (24) 

Denoting the average over the initial conditions by (•)p, one readily finds that the density 
(x)p is conserved: 

(x(x,t))p = ao + 2co. (25) 
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However, due to the Poissonian initial conditions which imply 

(a(x, 0)a(x', 0))p = + ao5^'^\x - x), (26) 

(c(x, 0)c(x', 0))p = cl + Co^^'^Hx - x), (27) 

one finds 

(X(x, t)2)p = (ao + 2co)' + (ao + 4co) (STrDt)-'^/^^ (28) 

This approximation gives satisfactory results concerning the power law approach to 
equilibrium, however it is unable to predict the correct sign of the correlations. As a 
consequence of this erroneous sign, the density of the C particles would always approach 
its stationary value from above, leading to a non-monotonic behavior, when cq < Cqo- 



2. 3. Approximation scheme for the concentration of C particles 

In this subsection, we would like to compute the approach of the density to its stationary 
value. Let us define &(x, t) = c(x, t) — c^o and 5\;(x, t) = x(x, t) — (ao + 2co). The 
Langevin equation for 5c is then given by 

(9t-DV' + a-AA)&(x,t) 

= 4A &(x, tf - 4A (5>^(x, t) &(x, t) + \ (^(x, tf + ^oaa - /i) (5%(x, t) (29) 

where we put oaa = 4Aaoo + /i- The explicit solution of (29) is unknown. However, we 
can obtain the large time behavior of (&), by exploiting that is a Gaussian random 
variable with vanishing variance when t goes to infinity. The formal solution of (29) can 
be written in the form 

6c = g^[6c]+g[A\6c^ -A\5x5c + X5x^ + \{aAp,- ^i)5xl (30) 
where to simplify our notation we have introduced 

g[f]M = |Jdt|dVe-'^-^(*-*')Go(x-x',t-t')/(x',t') (31) 

and 

go[f]i^,t) = J dVe-'^^*G'o(x-x',t)/(x',0). (32) 

The first term in (^) comes from the initial condition. It simply reduces to (cq — 
Coo)e-"^*. 

Iterating this solution will give eventually a series in ^o[&]) ^[^] and with 
appropriate insertions of the operator g. Three kinds of terms then occur: terms 
containing only power of ^olifc]) terms containing only power of g[Sx] and and 
mixed terms. One easily sees that the first kind of terms gives exponential decay in 
the long time limit, they can thus be discarded as we are interested in the asymptotic 
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time regime of (&) and After averaging, and due to the particular structure of 

the operator Qq, the mixed terms will also give exponential decay when t ^ oo. They 
can thus be discarded. Therefore, in the long time limit, (&) reads: 

(&) = \Qm')] - 2A(aAA - ^i)Q[{5xGM)\ + A(oM - fifG[{QWf)\ +■■■ (33) 

where the ellipsis stands for terms containing at least fourth power of 5x- In the following, 
we shall show that they give sub-leading contributions to the asymptotic time behavior 
of {5c). Let us now analyze the first term of that equation. It reads 

^[(<5^2^] = ^*dt'y'd'^xe-'^'^(*-*')Go(x-x',t-t')(5Y(x',t')') (34) 

The large time behavior of this expression can be obtained in several ways. One would 
be to integrate over the space dependence, and then to use the property 

*dt'e-"(*-*')/(t') = -/(t) + 0(/'(t)), (t-oo) (35) 
a 

which of course is only valid if f'{t) is negligible with respect to f{t) (in particular this 
property is false when f{t) is an exponential). In our case, one could safely use it, as 
(5\;2) (which plays the role of f{t)) decays with a power law. Thus one finds 

Qm')] = (36) 

Another equivalent way to obtain this result is to note that the leading behavior of the 
time integral is obtained when t' t, i.e. one can simply replace e"'^-^*^*"*'-' with the 
delta function 6[aAA{t — t')] and then use that 

G'o(x-x',0) =5('^)(x-x'). (37) 

In other words, in the long time limit, the operator Q can simply be replaced by a^. 
Applying this method for all the other terms of (33), one readily obtains the large time 
behavior for (&): 

(&) = f A - + ^(^M^) (^2) + . . . . (38) 



Note that the same result would have been obtained using the property ([351). The 
terms containing higher powers of Q[6x\ can be treated in the same way. They will give 
(subleading) contributions of order (Sx^)^, with n >2. Finally, one finds 

(&) = ^{6x') + ■■■ = ^(co - c^){87TDt)-'/' + ■■■ (39) 

"^AA '^AA 

Sub-leading corrections to that result are of order or t^^/^"^, and will generally 
depend on microscopic details such as the lattice constant, etc. 

In summary, using the fact that is a Gaussian variable with a vanishing variance, 
we have written & as a power series in 6x, in a systematic way (formally, using this 
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method, one can compute the sub-leading corrections, however the calculation may 
become rather tricky). The same method may be used to obtain (&^), yielding 

2 



2oaa 
oaa 



OAA 



(co-Coo)(87rZ}t) 



-d/2 



+ 



(40) 



Equation (40) implies that (&) is not a self-averaging quantity. Indeed, the relative 
density fluctuations are given by 



1/2 



Oaa 



U 2 I 1 — 1/2/ 

-^aa|co - Cool (»7r 



Dt) 



d/A 



(41) 



A/i2 

They diverge in the long time limit. As a consequence, it is very diflicult to conflrm 
numerically the validity of our result (39). Indeed, denote by (Ac) the average density 
measured in a typical simulation: 



(Ac) 



1 



2V 

Oaa 



[co-Coo)(87rZ}t)-^/2 



(42) 



(where L is the system size) 



(AcO 



\^2d 
1 

21^ 



Its fluctuations are 
d'^xjd'^x' &(x,t)&(x',t) 



^ jd'x (&(0,t)&(x,t)) 



Oaa 



Oaa 



(cq — Co 



so that the relative fluctuations reads 



(Ac^) - (Ac)' 



(Ac)^ 



1/2 



(oaa - 

2A/i2 



O^A I O-c 



871 Dt/L) 



d/2 



(43) 



(44) 



The fluctuations are negligible only if t <C L (in contrast with the diffusive behavior for 
which t <C I/^), a regime which is very difficult to obtain numerically. 

Intriguingly, the same results for the density and its correlations could be obtained 
by imposing {dt — DV^) & = in equation (29) and then solving the quadratic equation 
for §c. However this approximation is uncontrolled, and only the lowest order terms 
can be obtained. It is possible to explain why such a crude approximation works, by 
noting that even though & is a random variable, it should not vary too rapidly, as it 
only depends on 5x and not directly on the noise C,. 

Remark that the sign of (&) is given by {5x^) (or equivalently by cq — Coo)- This 
means that if initially {5c) is positive, so will it be for large time. No information is 
given for intermediate times, but a non-monotonic behavior would be surprising. 

The computation of the density correlation functions goes along the same lines. 
Using that a = a^o + 5x — '25c and equation (H), one obtains (when t oo) 

,2 



Ca{^,t) = al + 5^''\^) + -^(aoo - ao)(87rZ}t)-^/2/ -.V8i?* _ 2\5^''\^] 

Oaa ^ 



(45) 
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where the Gaussian factor e ^^/^-D* comes from the expression for (5\;(x, t)6x{0, t)) (which 
is easily obtained from (]T8|)). The correlation length C,a is given by 

/ d'^a;a;^[C*a(x, t) — a^] AdirDtlaoo — ao\ 



Sa 



^ 1 

ooj 



(46) 



The absolute values ensure the positiveness of ^,1 (note that Ooo — ao + ^ooCaa//^^ > 0). 
When ao > aoo, the second moment of Ca(x, t) — is negative, indicating that the A 
particles are negatively correlated. The same conclusion holds for the C particles. 



3. The two-species reversible reaction A + B ^ C 
3.1. The model 

In this section we study the reversible A + B ^ C reaction using the previous approach. 
As before, our starting point is the continuous time master equation describing the 
process on a c?- dimensional hypercubic lattice, allowing multiple occupancy. Let 
P({/}, {m}, {n}; t) be the probability to find the configuration {Z}, {m}, {n} at time 
t. The set {1} describe the occupation numbers on each lattice site for the A particles, 
{m} is used for the B particles and {n} for the C's. The master equation is then 

|p({/},{m},M;t) 

= ^ E E[(^e. + 1)P(. 1, /e. + 1, ... , W, {n}- 1) 
-kPi{l},{m},{n};t)] 

+ ^ E E[K, + 1)P({/}, . . . , m, - 1, m,, + 1, . . . , {n}; t) 
-miP({/},{m},M;t)] 

+ ^ E EtK + 1)^({0, {^}. . . . , - 1, + 1, . . . ; t) 

^ i ei 

-n,P({/},{m},M;t)] 

+ E[(^i + 1)K + . . , /i + 1, . . . , . . . , mi + 1, . . . , {n]- 1) 

i 

-mi{mi-l)P{{l},{m},{n]-t)\ 

+ E[(^^ + W, . . . , + 2, . . . ; t) 

-n,P({/},{m},M;t)]. (47) 



This master equation has the same structure than equation (4) , and the same notation 
is used {D}y is the diffusion constant of B particles). For the time being, we choose 
homogeneous initial conditions given by an uncorrelated Poissonian distribution on each 
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site and for each species: 



P({/}, {m}, M; 0) = e-«°-^°-<^° 11 77^^' (48) 

where 6o is the initial average occupation number of B particles. 
In the field theory formalism, the action is given by 

j d'^xj^* dt [a{dt - DaV^)a + b{dt - DhV^)b + c{dt - D^W^)c 

+ (ab + a + b-c) {Xab - /jlc) - 6{t) (ooo + bob + cqc)] , (49) 

where bo — bo/i^. In term of Langevin equations we find 

{dt - i^aV2)a(x, t) = -Aa(x, t)6(x, t) + ^c{^, t) + Ca(x, t) (50) 
{dt - DbV^)b{x, t) = -Aa(x, t)6(x, t) + ^c(x, t) + (^(x, t) (51) 
{dt - DcV^)c{:si, t) = Aa(x, t)6(x, t) - iic{x, t), (52) 

where and ^6 are two complex Gaussian noises with zero mean value and whose 
correlation are given by 

(Ca(x, t)Ca(x', t')) = (a(x, i)C6(x', t')) = (53) 

(Ca(x, i)C6(x', i')) = 2(/ic(x, i) - Aa(x, i)6(x, t))S^''\^ - x')S{t - t') 

= -2dt{c{t))5^'^\x - x')5(t - t') (54) 

(for homogeneous initial conditions). As for the A + A^C reaction, the noise vanishes 
at equilibrium. The equilibrium densities are thus given by their mean-field solution: 

(55) 

which together with the two conservation laws 

aoo + boo + 2coo = oo + 6o + 2co (56) 
Ooo - ^'oo = ao - ^'0, (57) 
permits us to obtain the equilibrium densities as a function of the initial conditions. 



3.2. Conserved quantities 

The approach toward equilibrium will be obtained through the same steps as before. 
First we write the Langevin equation for the conserved quantities and then, once these 
equations solved, we can set up a systematic approximation scheme for Sc. 

Let us introduce ip — a — b and % = a + 6-|-2c. Prom equation (50)-(52), one readily 
obtains the two following Langevin equations: 

(9t-L>VXx,i) = C^(x,t), (58) 
{dt-DV')x{^,t)^CxM, (59) 
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where in order to simplify, we supposed Da = = = D (the case of different 
diffusion constants will be considered in section The noises and have a vanishing 
mean and their two-point correlations are 

(C^(x,t)C^(x',t')) = -(Cx(x,t)Cx(x',t')) 

= 2at(c(t))5(^)(x - x')5(t - t'). (60) 

(g(x,t)C^(x',t')) =0 (61) 

The solution of these two Langevin equations has the same form as equation (^8|), with 
appropriate initial conditions (V'(x, 0) = — &o and x(x, 0) = Oq + &o + 2co). ^Z' and x 
have then a Gaussian distribution with a non vanishing mean 

(^) = ao - &o, (x) = ao + &o + 2co, (62) 
and their variance is 

ii^') - = -((X') - ix)') = 2 r dt,[8nDit - t.r'/'dMti))- (63) 

JO 

In the long time limit, we find again the power law decay: 

(^2) - iij)' = 2{c^ - co)(87^D^)-'^/^ {t ^ oo). (64) 

3.3. Approximation scheme 

The next step is to set up our approximation scheme for {&) = c — Coo- Let us define 
dip = ip — {aQ — bo) and Sx = X ~ {C'O + + 2co). The Langevin equation for & reads 

= A &(x, ty - A 6x{^, t) &(x, t) + \XM-^, tf - &/'(x, tf\ 

+ |(oab - /u) (5\:(x, t) - \\{ao - bo) &/'(x, t), (65) 

where we put oab = A(ao + &o + 2co — 2coo) + This equation possesses the same 
structure as equation (29). By repeating the same scheme, we may then obtain the 
large time behavior of (&). Finally, we find 

ao - bo^ 21 



1- A 



40-AB 



((V ) + ■ ■ • (66) 



or, by putting the large time expression for {Sx ) and {Sip ) 

(&) = l^WL + - A'(«o - &o)'](co - Coo)(87rDt)-^/2 + • • • (67) 

Note that + /i^ — A2(ao — &o)^ = 2/i(oAB + 2Acoo) > 0. 
Similarly, one obtains the C-particle correlation: 

(^e) = ^^(co - c^)(87rDt)-'^/2 + • • • (68) 
Sub-leading corrections to these laws are of order t^'^ or t^^^"^^^. 
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Extensions 



4.I. Correlations in the initial condition 



Recently, Yang et al showed that in the A + B ^ C case, and for correlated 
initial conditions, the power law approach ( ]68| ) could be modified. More precisely, for a 
particular initial condition, they found that in one and three dimensions the approach to 
equilibrium was of the order i.e. faster than in the uncorrelated case. Whereas 

it seems to be at odds with the fact that the power law should be universal (see the final 
section), we shall show in this subsection how such a behavior can be obtained within 
our approach. 

Let us consider the reaction A + A ^ C , when initially a fraction of the total density 
of the A particles are not distributed independently, but are disposed in pairs having a 
separation radius a (each A- A pairs are supposed to be independently distributed). We 
still denote the total density by ao? and the density of pairs will be denoted by uq. In 
order to take this condition into account in our formalism, one should add the following 
term in the action (7) 

no 



j d-^x j d^y a(x, 0)a(y, 0)5(|x - y| - a). (69) 



(srf is the surface of d dimensional sphere of radius unity). In the language of Langevin 
equations, this new term translates into a new contribution to the noise correlation, 
namely 

(C(x, t)C(x', t')) = -29(c(t))5W (x - x')5(t - t') + 25(|x - x'| - a)6{t)6{t'). (70) 

The only relevant effect of this new contribution is to affect the variance of Xi to which 
one should add 

{^')corr = j^'^^'j dVGo(x-x',t)Go(x-x",t)5(|x-x'| - a) 

= 2noexp{-3a^/ADt){8nDt)-'^/\ (71) 
For long time, {Sx^) becomes then (we used that 2(co — Coo) = ctoo — cto) 

{6x') = K - (ao - 2no)](87rDt)-'^/2 ^ . . . (72) 

oq — 2no represents the initial density of uncorrelated A particles. Inserting this result 
into formula (39) one obtains the C-particle density: 

(&) = -^(aoo - ao + 2no)(87rZ}t)-'^/2 ^ . . . (73) 

Initial correlations could then lead to non-monotonic behavior of (&) , as aoo — clq and 
Ooo — + 2no do not have necessarily the same sign. 
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A case of special interest occurs when we choose 2no to be exactly a^o — cto- 
The amplitude of the leading term of the approach to equilibrium vanishes and one 
should consider the sub-leading corrections to {5x^)- As a consequence the approach 
to equilibrium will be faster. It is easy to verify that these corrections will be of order 
j.-d/2-i ^ with an amplitude that depends on microscopic parameters such as a, the initial 
correlation length, or £, the lattice constant. Note that the condition 2no = Ooo — do can 
be easily obtained experimentally: consider a system where A and C particles are at 
equilibrium. At time t = 0, we excite some fraction of C particles such that they break 
up into pairs of correlated A particles (this excitation can be for example obtained by a 



photo-flash [^). Let us denote by and Cqo the concentrations of A and C particles 
before the excitation (note that we were at equilibrium: Aa^ = /iCoo)- At time t = 0, the 
initial densities will be = + 2no and Cq = Cqo — n^. Now, as Oq + 2co = cioo + 2coo, 
one must then have a^o = floo, which implies 2nQ = — Ooo- This initial condition 
correspond exactly to the case studied by Yang et al pO|, for the two-species reaction. 

Similar conclusions can be drawn for the A + B ^ C reaction when initially a given 
amount of A and B particles are distributed in pairs of radius a (initial correlations 
among particles of the same species will not be considered, as they should decay 
exponentially fast). The action get modified in the same way as previously, by adding 
a new initial term: 
no 



j d^x J d% a(x, 0)6(y, 0)5(|x - y| - a). (74) 

In terms of the Langevin equation it also modifies the noise (, leading to a new 
contribution for both {Sip"^) and {Sx^) (see (71)). Finally, one gets 

{6x') = = K - («o - %)](87rZ}t)-'^/2 ^ . . . (75^ 

Here uq is the initial concentration of correlated pairs, it is also the initial density of 
correlated A (or B) particles. Contrarily to the irreversible A+B case, the presence 
of correlations in the initial state does not modify the long time behavior (except in some 
special limits as discussed above). 



4-2. Unequal diffusion coefficient 

The case of unequal diffusion constant is of special interest, because if, for some reasons, 
one kind of particles moves very slowly compares to the other, one could legitimately 
question the validity of our previous results. In fact we shall see that even in the worse 
case (one species at rest) the power law does not change, only the amplitude is affected. 

Let us consider first the A + A ^ C reaction, the generalization to the A + B ^ C 
reaction being straightforward. The two Langevin equations associated with that model 
are 

{dt - DaV) 6x + 2{Da - D,)V' & = C, (76) 



17 

{dt - D,V^ + o-aa) (fc = 4A(fc2 _ 4A5^&+ A(5%2 + l{a^ - ^) (77) 

The first consequence of liaving two different diffusion constants lies in tlie fact tliat tlie 
Langevin equation for Sx is no longer closed, but contains a term proportionnal to Sc. 
The main question we want to address, is how does this extra term affect the large time 
behavior of the One quick and false answer would be to say that, as & is a slowly 

varying random variable for large time, it should not modify the large time behavior 
of Sx- This picture is however not true, as & itself depends on in a non trivial way. 
In fact, the better way to treat this term is to insert it into the propagator, which is 
then no longer diagonal. Putting also in it the term proportionnal to 6x in equation 
(77) (which was previously treated like an interaction term), the propagator is given, in 
Fourier and Laplace transform representation, by the inverse of the following matrix 

V-2(^AA-/i) S + DcP' + CTaaJ 

where p and s are the Fourier and Laplace conjugated variables of space and time (this 
last expression is very easily obtained in the field theory formalism by considering the 
quadratic terms). For the Sx^x propagator, one then finds in the (p, s) representation 

^ r ^ g + DcP'^ + q-AA 

The 6xSc propagator is given by 

G^MlP, - ^ j^^p, ^ ^^^^^ ^ ^^^2) + (^^ _ jj^)^p2- 

The two other propagators could be obtained in the same way, but we are not interested 
by them. Our purpose now is to compute the large time limit of (Sx^)- Once this 
obtained, we shall then use our approximation scheme (applied to equation (77)), in 
order to derive its large time behavior. 

For computing the second moment of 6x, it is easier to consider a diagrammatic 
representation (in this paragraph, we follow the line of reasoning developped in for the 
irreversible reaction A + B 0, when Da 7^ -De)- A dx^x propagator will be represented 
by a dashed line and a && propagator by a full line. Off-diagonal propagators will be 
represented by mixing of the two lines (see figure 0). From the Langevin equation (77), 
three different vertices (four with the noise) can be identified, their representation is 
also given in figure |^. As we shall eventually average over the noise, only diagrams 
containing two merging noise lines at their beginning can subsist. To obtain the second 
moment of 6x, we simply have to draw all diagrams ending with two merging 6x lines (see 
figure If Da = Dc, only the first diagram would give a non- vanishing contribution (in 
agreement with equation (|20|)). In fact, the subsequent terms of the sum all contain at 
least one of the three sub-diagrams shown in figure ^ However, in the language of field 
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time 




Figure 1. Set of propagators (upper part) and vertices (lower part) needed for the 
diagrammatic representation of the Langevin equations. The noise term is represented 
beginning with a cross. The arrow on the left represent the direction of the time. 

m= C['y + <^:~-, + ^'Cil +■•• 

Figure 2. Diagrammatic representation of (Sx^). The dot stands for noise- noise 
correlation. 

,«v «^ 

Figure 3. Effectives vertices appearing in all diagrams besides the first one of figure ^. 
They all give sub- leading contribution to {Sx^}- 



theory, these sub-diagrams give rise to effective vertices of the form SxV'^Sx'^, 5\;V^5\;(fc 
and 5\;V^&^, which, by simple power counting, turn out to be irrelevant. They will only 
give sub-leading contribution to (Sx^)- Hence, its leading term is only given by the first 
diagram of figure 

To compute this contribution one first needs to obtain Gs^SxA (P? ^) space. 
Inverting the Laplace transform, one finds 

X {[OAA - (Da - D,)p^] sinh[ir(p)t] + r(p) cosher {p)t]} 6 (t) (81) 

where 

r{p) = ((o-AA + {Da - D,)p^f - Afi{Da - D,)p^y^\ (82) 
One is now in position to obtain (Sx^)- To leading order, it is given by 

{Sx') = -2 /Jdt'/^[G^^,o(p,t-Of5tW)), (83) 
whose large time behavior (which can be obtained by taking the large cm and /i limit) 
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reads 

{dx') = -2{8nD,s)-'^' fdt' (t - tr"'dt{5c{t')) + • • • (84) 

with 

D,s = D, + {Da-D,)^ (85) 

This is nothing but the result of equation (pO|), with D replaced by D^s- 

We can now come back to the Langevin equation for 5c which is exactly the same 
than the one written for equal diffusion constant (equation (29)), with D replaced by 
Dci and then use our approximation scheme as before. The final expression for (&) will 
then be given by equation (39) with D replaced by D^q. 

It is particularly interesting to consider the cases when either Da or D^. vanish. If 
Dc = 0, one finds D^s = DafJ^/ojs^- Although the C particles do not move, one still 
observes a power law decay of the concentration, but with a smaller diffusion constant 
than when Dc = Da {fi < oaa)- The fact that the C particles are at rest is compensated 
by the movement of the A particles, which leads to an effective diffusion of the C 
particles. More surprising is the case Da = 0. At first sight, one could expect that 
because the A particles do not move, the forward reaction is essentially inoperative 
(only A particles at the same site could react), however, one should not forget that the 
C particles still move and that they effectively carry two A particles (allowing thus some 
mixing of the A particles), and secondly this motion still allows the fluctuations of 6x 
to be smoothened diffusively. 

For the A + B ^ C reaction, the case where Dc differs from Da = Db can be studied 
in exactly the same way. All we said for the A + A ^ C reaction is still valid, with 
the modification that, as 6ip only depends on the random variables a and b, it is not 
modified. Only 6x changes. As a consequence, the formula (p^ ) is still valid, and for 6x 
one finds 

{5x^) = 2(coo - Co)(87rDefft)"'/', (t ^ oo), (86) 

where Des is given by equation (|85D with oaa replaced by ctab- The expression for (&) 
becomes then 
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2 / , (Tab- fiDa- Dc\ '^1'^ ^ 



V om Da ) 
xico-CooKSnDaty/^. (87) 



When Dc = 0, the amplitude is slightly modified (compared to the case Dc = Da) but 
one could still consider the diffusion constant to be Da- Note that the case Da = 
cannot be treated by this formalism. 

The case where all the diffusion constants are different can be considered as well. 
We shall not treat it here, but the result should not differ too much from the previous 
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one. Indeed the main new ingredient is that {Sip"^) is modified. However it has been 
shown that, for the A + i? — reaction, this leads only to a change in the amplitude 
P]. We expect this result to be only slightly modified in the reversible case. For dx 
the same analysis as before applies, but with more complicated expressions. The cases 
Da = or Dfe = cannot be treated with this formalism. 



4-3. Segregated initial conditions 

Our approach can be extended as well to other initial conditions. The A + B ^ C 
reaction with initially segregated reactants (say, the A particles on the right, the B ones 
on the left and no C particles) is of particular interest, mainly due to the dynamics of 
the reaction front. In the irreversible case, it is now well established that the width 



of the front increase with a power law w{t) ~ |^T], |2^, Two different cases 

may be distinguished: above two dimensions, the exponent takes its mean-field value 
a = 1/6, whereas below two dimensions fiuctuation effects play a dominant role, leading 
to a = l/[2{d + 1)]. Extrapolating our previous results, one could expect that the 
fluctuations in the conserved quantities will play an important role. In fact it appears 
that the width of the front increases indeed with a power law w{t) ~ t", but this 
exponent takes its mean-fleld value (a = 1/2) for any dimensions. Thess surprising 



results have already been obtained by Chopard et al using scahng arguments and 
numerical simulations. In this subsection, we shall show how this behavior can be 
conflrmed within our formalism. 

The problem is described by the Langevin equations (50)-(52) with the initial 
conditions 

V>(x,0)=no[^^(a;i)-e(-Xi)] (88) 
x(x,0)=no (89) 
c(x,0) = 0. (90) 

For simplicity we have chosen the particle to have the same diffusion constants; moreover 
both reacting species A and B are supposed to be homogeneously distributed (with 
density uq) in their respective semi-inflnite sub-space. Let us denote by c the mean-fleld 
C particles density. One easily shows that in the long time limit, it takes a scaling 
behavior: 

c(x,t) ~Coo(0, (91) 



where ^ = xi/yADt. The exact form of Coo(0 is unimportant. It can be obtained by 
equating to zero the right hand side of the rate equation for the C particles. Deflning the 
width of the front by the square root of the second moment the density of C particles, 
one immediately obtains, in the mean-fleld case, the exponent a = 1/2 for the reaction 
front. 
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In order to take into account the fluctuations, one needs first to integrate the 
equation for the conserved quantities and x- the long time limit, one can show 
that 

(<502) = _(^sx') ^ -{STiDt)-"/' r^^e-^\^{^ - y). (92) 

J-oo y TT 

The Langevin equation for the C particles can be rewritten using & = c — c, so that 

[dt - L»V']& = \X{5x^ - (50') + - X5x5c 

+ \[(7{xi,t) - ii]5x. - |Anoerf(^) + (T{x-i,t)5c (93) 

where (T{xi,t) = Xhq — 2Xc{xi,t) + /x, and erf is the error function. This equation is 
very similar to equation (65), with the difference that a has become time and position 
dependent. As a consequence, we are unable to write in a closed form the propagator 
related to this equation, and thus to apply our approximation scheme. However we can 
still deduce the large time behavior of the density (&) by equating the right hand side 
to and solving the quadratic equation in This crude approximation has proven to 
give accurate results in the homogeneous one- and two-species reactions. In addition, 
the same results arc obtained by assuming (on physical grounds and by examining its 
differential equation) that the propagator behaves like 

G{x,x',t,t') ^exp[-a{xi,t){t-t')]Go{x-x',t-t') (94) 

when t' ^ t (note that in the homogeneous case, this limit gave the asymptotic time 
behavior) . 

Finally, one finds 

(c) = c{x,,t) + — ^{/x' + ^(^1,^)' - X^nll^^mfm')- (95) 

As in the homogeneous case, the mean-field asymptotic solution is approached with a 
power law. From the last equation one immediately obtains that the width of the front 
is governed by its mean-field exponent a = 1/2. This result is easily explained: the 
spreading of the front is given by the diffusion of the C particles. Moreover, once the 
backward reaction is allowed, it has been shown that the C particles will always diffuse 
with a non- vanishing effective diffusion constant (even when Dc — Q). Hence the width 
of the front should grow like the square root of time, independently on the fluctuations 
which are governing only the approach to the equilibrium, and not the spreading of the 
C particles. 

5. Discussion and concluding remarks 

As we mentioned in the beginning, our model allows for multiple occupancy of each site 
and contains only single site reactions, a property which has considerably simplified our 
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analysis (leading in particular to the simple form of the Langevin equations). It is thus 
natural to question about the universality of our results. 

To answer this question, let us consider the following two Langevin equations 

{dt - DV)^ = C (96) 
{dt - DV)$ = (97) 

i+j>l 

with {0 = and 

(C(x,t)C(x',tO) =r(t)5('^Hx-x')5(t-t') (98) 

(homogeneous case). By writing these two equations, we have implicitely assumed that 
initially "^{x, 0) = and that limt^ooi^) = (otherwise a constant term a^Q ^ should 
be added). From equation (96), one immediately obtain in the long time limit 

(^2) = (87rD)-'^/2 /■* _ ijydli^ijJ^ (99) 
h 

(to cure the divergence when i! — > t, the integrand should be multiplied by a cut-off 
function), leading to the following long time behavior for ($) 

($) = -^(a^.aao - aoiaiiaio + ao2a?o)(87rD)-^/2 /*dt' (t - t')"'/'r(t'). (100) 
agi JO 

If aoi 7^ 0, the leading behavior is given by /q dt' (t — i!)^'^l'^V{t') (as long as 
0^01^20 ~ c^oi'^ii'^io + '3o20io 7^ 0). If /q°° dt' F (t') is finite (this implies in particular 
that the noise dies out at equilibrium), one recovers the t^'^l'^ power law. This shows 
us that in general the power law does not depend on the structure of the equation for 
but only on the presence of \Ef, i.e. of a diffusive conserved mode in the model. The 
condition aoi 7^ (here aoi plays the role of oaa or oab) implies the presence of a non- 
vanishing mass in the field theory for the field $. This is responsible for the presence of 
an exponential decay in time of the propagator, which, in turn, plays a central role 
in the large time behavior. 

Equation (|100|) shows that the amplitude depends only on the most relevant 



operators (the value of the coefficient ajj, with 1 < z + j < 2). In particular, a model 
with an exclusion principle on each lattice site can be handled in the same way. Whereas 
the exclusion condition will give rise to new contributions for the noise correlation and 
for various operators, one can reasonably expect on physical grounds that the conserved 
quantity will behave diffusively, leading to the t~'^l'^ power law. However, the equation 
for the evolution of $ will be modified, leading to new expressions for the steady-state 
densities (which are again simply obtained by solving the corresponding mean-field 
equation) and the amplitude. Note that the exact expression of the equilibrium densities 
depends explicitly on the complete equation for (the unshifted version of) $, but as the 
noise dies out at equilibrium, it is still given by a mean-field equation. 
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The addition of a new diffusive conserved mode (like in tfie case of the A + B ^ C 
reaction) clearly does not alter the power law, but the amplitude. A generalization of 
our analysis is also straightforward for the reversible aggregation + An Am+n, 
for which the total number of monomers A is conserved. It is however not possible 
to extend our results for the reversible coagulation process A + A ^ A, for which no 
such conservation law occurs. In fact this last reaction can be described by an action 
equivalent to the one obtained for directed percolation, in the active phase. One expects 
thus the upper critical dimension to be 4. 
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